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Black hole accretion disks in the canonical low-hard state 
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ABSTRACT 

Stellar-mass black holes in the low-hard state may hold clues to jet formation and ba- 
sic accretion disk physics, but the nature of the accretion flow remains uncertain. A stan- 
dard thin disk can extend close to the innermost stable circular orbit, but the inner disk may 
evaporate when the mass accretion rate is reduced. Blackbody-like continuum emission and 
dynamically-broadened iron emission lines provide independent means of probing the radial 
extent of the inner disk. Here, we present an X-ray study of eight black holes in the low-hard 
state. A thermal disk continuum with a colour temperature consistent with L oc T A is clearly 
detected in all eight sources, down to «5xl0 -4 LEdd- In six sources, disk models exclude a 
truncation radius larger than 10r g . Iron-Ka fluorescence line emission is observed in half of 
the sample, down to luminosities of r*1.5x 10 -3 LEdd- Detailed fits to the line profiles exclude 
a truncated disk in each case. If strong evidence of truncation is defined as (1) a non-detection 
of a broad iron line, and (2) an inner disk temperature much cooler than expected from the 
L oc T 4 relation, none of the spectra in this sample offer strong evidence of disk truncation. 
This suggests that the inner disk may evaporate at or below « 1 .5 x 10 -3 LEdd- 
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1 INTRODUCTION 

The innermost regions of accreting black holes and the underlying 
physics of the accretion process can be successfully studied using 
the X-ray spectra of black hole binaries. Over the past years, such 
studies have revealed that Galactic black hole binaries (GBHBs) 
radiate in various distinct spectral states characterised by the rela- 
tive strength of their soft and hard X-ray emission (McClintock & 
Remillard 2006). 

The X-ray spectrum of black hole binary systems in their low- 
hard state is dominated by a powerlaw with relatively low luminos- 
ity (^0.05 LEdd', Maccarone 2003), a photon index Y in the range 
~1.4-2 and an exponential cut-off at about 100 keV. This is in con- 
trast with the high-soft state, where the spectrum is dominated by 
a quasi-black body component with a characteristic temperature of 
~1 keV, and the highly luminous very-high state (L~LEdd) with 
a powerlaw component having a flux comparable to that of the soft 
blackbody. In the latter cases, the powerlaw does not show evi- 
dence for a high-energy roll-over. Superimposed on the continuum 
is the presence of various reflection features (Ross & Fabian 1993; 
Miller 2007) with the dominant component being the Fe-Ka emis- 
sion line. In the inner regions of the accretion flow the reflection 
spectrum appears blurred due to the combination of various rela- 
tivistic effects. The degree of gravitational blurring is strongly de- 
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pendent on the inner radius of the reflecting material (Fabian et al. 
1989; Laor 1991). 

Theoretically, the soft, quasi-black body component observed 
in the high-soft and very-high state is generally agreed to origi- 
nate in a standard, optically thick, geometrically thin accretion disk 
(Shakura & Sunyaev 1973) extending to the innermost stable cir- 
cular orbit (ISCO). For a non-rotating, Schwarzschild, black hole 
this radius is equal to 6r g where r g = GM/c 2 , and decreases to 
«1.2r g for a maximally-rotating, Kerr black hole (Thorne 1974). 
The powerlaw component dominating the X-ray spectra in the low- 
hard state, on the other hand, is believed to be produced by the 
inverse Compton scattering of soft photons in a thermal, optically 
thin region (Shapiro et al. 1976; Sunyaev & Titarchuk 1980). How- 
ever, there is no general consensus on the geometry of this optical- 
thin region. In the accretion disk-corona model the hard powerlaw 
comes from either a patchy corona, possibly powered by magnetic 
flares (Di Matteo, Celotti & Fabian 1999; Beloborodov 1999; Mer- 
loni, Di Matteo & Fabian 2000; Merloni & Fabian 2001) or the 
base of a centrally located jet (Merloni & Fabian 2002; Markoff 
& Nowak 2004; Markoff, Nowak & Wilms 2005). In both cases 
the thin accretion disk extends close to the ISCO. An alternative 
model has the thin disk truncated at large distances from the black 
hole (Narayan & Yi 1995; Esin et al. 1997) with the central region 
replaced by an advection-dominated accretion flow (ADAF). 

In both interpretations for the geometry of the accretion disk 
in the low-hard state it is expected that the thermal-disk component 
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will have a low effective temperature. The low mass accretion rates, 
and corresponding low luminosity, observed in systems in the low- 
hard state implies a peak disk temperature < 0.4 keV. In the case of 
a geometrically thin disk extending close to the ISCO, this tempera- 
ture closely follows the L <x T 4 relation. The expected temperature 
departs dramatically from this relation when the disk is truncated at 
the distances predicted by the ADAF models (see e.g. McClintock 
et al. 2001). A thermal component has been observed in a number 
of sources in the low-hard state. In particular, XMM-NewtonlRXTE 
observations of GX 339-4 (Miller et al. 2006; Reis et al. 2008) and 
Swift J1753.5-0.127 (Miller, Homan & Miniutti 2006; Reis et al. 
2009b) as well as Chandra observation of XTE Jl 118+480 (Reis, 
Miller & Fabian 2009) have shown that the accretion disk in these 
sources are consistent with extending close to the ISCO. Evidence 
for the disappearance of this component has usually been the re- 
sult of observations made with RXTE which lacks the low energy 
coverage. 

As mentioned above, a further feature in the X-ray spectra of 
BHBs are the various reflection signatures associated with the re- 
processing of hard radiation by the cool accretion disk. The most 
prominent of these features is often the broad, skewed Fe-Ka line 
observed in a number of BHBs (Miller 2007, 2009), Seyfert galax- 
ies (Tanaka et al. 1995; Fabian et al. 2009) and accreting neutron 
stars (Cackett et al. 2008; Reis, Fabian & Young 2009). New fits 
to some ultra-luminous X-ray sources (ULXs; Caballero-Garcia & 
Fabian 2009) spectra suggest that a disk reflection may also be ob- 
served, though more sensitive spectra are needed to confirm this 
possibility. The strength of the reflection features can be typified by 
the equivalent width of the iron-Ka line, Wkc which is expected 
to correlate linearly with the reflection fraction R = Q,/2n, where 
Q. is defined as the solid angle covered by the accretion disk as 
viewed from the hard X-ray source. George & Fabian (1991) have 
shown that for a neutral accretion disk with solar abundances the 
reflection fraction closely follows R « Wkc/ '180 eV. For stellar 
mass black holes in the low-hard state, R is observed to be below 
~ 1. This "weak" reflection can be interpreted as either a reces- 
sion in the accretion disk, as in the ADAF interpretation (Esin et 
al. 1997), a highly ionised inner disk (Ross, Fabian & Young 1999) 
or mildly relativistic motion of the hot corona away from the disk 
(Beloborodov 1999). Using the latter model, Beloborodov showed 
that reflection fractions as low as ~ 0.3 can be be obtained in the 
low-hard state without invoking disk truncation. 

In the following chapters we will present a systematic analysis 
of a sample of eight stellar mass black hole binaries in the low-hard 
state with the goal of better understanding the accretion disk and 
flow geometry of such systems. We find that for the sources ob- 
served in this work the accretion disks are consistent with extending 
close to the ISCO and suggest a set of observational criteria for the 
support of disk truncation. We start in §[2] with an introduction to 
the various systems investigated. This is then followed by a detailed 
analyses of the property of the thermal component (§§[3TTJ|3T4j and 
the reflection component ( i|3,5l >. Our results are summarised in §|4] 



2 OBSERVATIONS AND DATA REDUCTION 

Our aim is to present a selection of high quality spectra of various 
X-ray binaries in the low-hard state. In such a state the temperature 
of the accretion disk is usually observed to be below ~0.4 keV and 
for this reason a strong requirement of our selection is that the data 
extends to a similar energy range. Observations made solely with 
the Rossi X-ray Timing Explores (RXTE) are thus not used in this 



work due to its low energy cut-off of approximately 3 keV. The 
bulk of the work presented here utilises observations made with 
either XMM-Newton or Suzaku, with two further sources observed 
with Chandra and SWIFT. In what follows we introduce the various 
sources which meet this criterion and describe the observations and 
data reduction in detail. 

2.1 GX 339-4 

GX 339-4 is a dynamically constrained black hole binary located at 
a distance of 8 ± 2 kpc (Zdziarski et al. 2004). Although its mass 
is not yet known, it is likely that GX 339-4 is amongst the more 
massive of the stellar mass black hole sources. The mass function 
of the system has been constrained to ~6 Mq (Hynes et al. 2003 
and Munoz-Darias et al. 2008) and results derived from radio ob- 
servations suggests a low (8 ^ 30°) inclination for the inner disk 
(Gallo et al. 2004). Similar low inclinations have also been found 
when modelling the reflection features present in the spectrum of 
GX 339-4 in outburst (Miller et al. 2004). The apparent low in- 
ner disk inclination mentioned above does not necessarily imply a 
massive (> 30 Mq) stellar mass black hole since it is possible that 
the inner disk is warped similarly to that observed in GRO J 1655- 
40 (Martin, Tout & Pringle 2008) and V4641 Sgr (Martin, Reis & 
Pringle 2008). For the purpose of this work we assume a black 
hole mass in the range of 10-20 Mq and an inner disk inclination 
between 10-30 degrees. We note here that throughout this paper 
we will adopt the largest range for the physical parameters (mass, 
distance and inclination) available in the literature for the various 
sources, unless the parameters have been well constrained. 

Prior analyses of GX 339-4 in both the very-high (Miller et al. 
2004) and intermediate (Miller et al. 2008) state have shown that 
the accretion disk extends close to 2r g in those states. A similar 
conclusion has been reached by jointly fitting the reflection fea- 
tures present in the XMM-Newton spectra of both the low-hard and 
the very-high state (Reis et al. 2008). In this work, the low-hard 
state spectrum of GX 339-4 described in Miller et al. (2004) is re- 
analysed using the more stringent reduction procedure detailed in 
Reis et al. (2008). We restrict our spectral analyses of the XMM- 
Newton EPIC-MOS data to the 0.5-10.0 keV energy range. The si- 
multaneous RXTE PCA spectrum is fit in the standard 3.0-25.0 keV 
energy range with 0.6 per cent systematic error and an edge at 
4.78 keV (r = 0.1) to account for a strong XeL absorption fea- 
ture. The HXTE spectrum is fit in the 25.0-100 keV range. When 
fitting the RXTE spectra along with that of XMM-Newton, a joint fit 
is achieved by allowing a normalisation constant to float between 
the various spectra. All parameters were tied between observations 
unless stated otherwise. XSPECv 12.5.0 (Arnaud 1996) was used 
to analyse all spectra presented in this paper. The quoted errors 
correspond to a 90 per cent confidence level for one parameter of 
interest unless stated otherwise. 

2.2 XTEJ1650-500 

Optical observations have constrained the mass function of 
XTE J1650-500 (hereafter J1650) to 2.73 ± 0.56 M (Orosz et al. 
2004). Together with a lower limit on the inclination of 50° ±3° 
this sets an upper limit to the mass of the central source of ap- 
proximately 7.3 Mq. A recent analysis made by Shaposhnikov & 
Titarchuk (2009) has placed a strong constraint on the mass of the 
central black hole of 9.7 ± 1.6 Mq based on the scaling prop- 
erties of quasi-periodic oscillations (QPOs). The discrepancy be- 
tween their results and that of Orosz et al. (2004) are attributed 
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to a possible mis-classification of the spectral class of the optical 
companion. Due to the uncertainty in the mass of J1650 we conser- 
vatively assume that it lies between 5.3 and 11.3 Mg. The distance 
to J1650 has recently been estimated at 2.6 ± 0.7 kpc (Homan et 
al. 2006). 

XMM -Newton observed J 1650 during the transition from the 
low-hard to the very-high state in 2001 (Miller et al. 2002). The 
source was mistakenly identified as being in the very-high state 
due to the presence of an unusually strong (at the time) iron line 
emission. However, the cold (~0.3keV) thermal component and 
low luminosity likely place the observation in a rising phase of the 
low-hard state (see e.g. Rossi et al. 2005). The presence of a broad, 
asymmetric Fe emission line led Miller et al. to suggest that the 
black hole in J1650 is rapidly rotating. This was later confirmed 
by Miniutti, Fabian & Miller (2004) using BeppoSAX observations 
taken both before and after that of XMM-Newton. Contrary to this 
interpretation, Done & Gierlinski (2006) showed that the bright, 
low-hard state, BeppoSAX observation is also consistent with a disk 
truncated at ~10-20r g if there is resonance iron-K absorption from 
an outflow disk wind with a velocity of ~0.15c. 

Here, we analyse the same XMM-Newton EPIC-pn observa- 
tion described in Miller et al. (2002). The source was observed in 
the "burst" mode with the "thin" optical filter in place. The spec- 
trum was extracted using a box region centred on the source with a 
20 pixel width. Background events were extracted from a 10 pixel 
wide box adjacent to the source. Single and double events were in- 
cluded in the analysis. Response files were created in the standard 
way using the tools rmfgen and arfgen. The FTOOL GRPPHAwas 
used to require at least 20 counts per bin in all XMM-Newton obser- 
vations presented in this paper. We restrict our study of the XMM- 
Newton EPIC-pn data to the 0.6-10.0 keV energy range. 



2.3 CygnusX-1 

CygnusX-1 (hereafter CygX-1) is amongst the better known and 
studied stellar mass black holes. Having been the first object gen- 
erally recognised as a black hole (Tananbaum et al. 1972) it has re- 
ceived considerable attention. The distance to the source was very 
early estimated at 2 kpc (Murdin & Webster 1971). This has been 
confirmed more recently by Massey et al. (1995) where a distance 
of 2.1 ± 0.1 kpc is found. The orbital inclination of the system 
is thought to be in the range of 25-50 degrees (Gierlinski et al. 
1999). The uncertainty in this value follows the fact that the ac- 
cretion disk in CygnusX-1 is thought to be precessing (Stirling et 
al. 2001; Romero et al. 2002) with an average angle to the line of 
sight of ~30° (Fender 2001). The mass of the black hole is also 
highly uncertain, partially due to the uncertainty in the inclination. 
Although it is generally assumed to be close to 10 Mg, its value 
has been found by various authors to range from 7-25 Mg (e.g. 
Ziolkowski 2005; Shaposhnikov & Titarchuk 2009). 

Young et al. (2001) showed that the various spectra of 
Cygnus X-l in both the low-hard and the high-soft state are consis- 
tent with the relativistic blurring of a thin, ionised accretion disk ex- 
tending close to the ISCO. This interpretation was later challenged 
by Barrio, Done & Nayakshin (2003). Using the high energy cover- 
age of RXTE the authors were unsuccessful in measuring the high 
energy break assumed to be characteristic of the ionised accretion 



disk modeQ and thus argued for a truncated disk interpretation of 
the spectra. 

Makishima et al. (2008) observed the low-hard state of 
Cygnus X-l in 2005 October with Suzaku. Although the data suf- 
fered heavily from both pile up and telemetry saturation, the au- 
thors found the presence of a strong soft component with a temper- 
ature of ~0. 19keV which was associated with an accretion disk 
extending to within ~250km. Assuming a mass of 10 M© this is 
approximately 16r g . We report on two recent Suzaku observations 
where the effect of pile-up is less prominent. 

Cygnus X- 1 was observed with Suzaku on two occasions on 
2009 April 8 (Obs ID 404075020; hereafter Obs 1) and April 14 
(Obs ID 404075030; hereafter Obs 2) for approximately 50 and 
33 ks respectively. The three operating detectors constituting the 
X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) were oper- 
ated in the "burst" clock mode in both observations and employed 
the "1/4 window" mode. The front illuminated detectors (FI; XIS0 
and XIS3) were operated in the 2x2 and 3x3 editing mode, whereas 
the back illuminated detector (BI; XIS1) was operated in the 3x3 
editing mode only. A total FI exposure of ~10.7 and ~8.6ks for 
Obs 1 and 2 respectively were available for each XIS camera. The 
corresponding BI exposures were 5.3 and 4.4 ks. Using the latest 
HEASOFT v6.6.2 software package we reprocessed the data from 
the Version 2 processing following the SUZAKU Data Reduction 
GuidH We started by creating new cleaned event files using the 
tool "xispi" and the script "xisrepro" as well as the associated 
screening criteria files. XSELECT was then used to extract spec- 
tral products. In both observations, source events were extracted 
from a region composed of the intersection of a rectangle of size 
150x1024 pixels with an annulus having 30 and 200 pixels inner 
and outer radius respectively. Background events were extracted 
from a circular region with a radius of 90 pixels away from the 
source. We used the script "xisresp'0 with the "medium" input to 
obtain individual ancillary response files (arfs) and redistribution 
matrix files (rmfs). "xisresp" calls the tools "xisrmfgen" and "xis- 
simarfgen". Finally, we combined the spectra and response files 
from the two front-illuminated instruments (XIS0 and XIS3) us- 
ing the FTOOL ADDASCASPEC. This procedure was repeated 
for each observation resulting in a total of four XIS spectra. The 
FTOOL GRPPHA was used to give at least 500 counts per spectral 
bin in every Suzaku observation mentioned in this paper. The Hard 
X-ray Detector (HXD) observations of CygX-1 will not be used in 
this work. The FI-XIS spectra were fit in the 0.7-10.0 keV energy 
range. The results obtained with the BI spectra are not presented in 
this work as it is not expected to be different to that of the FI. In 
all fits to Suzaku data presented in this paper the 1.9-2.5 keV en- 
ergy range was ignored due to the possible presence of unmodelled 
instrumental features. 

2.4 SWIFT J1753.5-0127 

Using multiwavelength observations, Cadolle-Bel et al. (2007) esti- 
mated the distance to SWIFT J1753.5-0127 (hereafter J1753) to be 
in the range of 4-8 kpc. This was based on the derived hydrogen 

1 If the accretion disk is ionised solely due to the hard X-rays from mag- 
netic flares then a high energy break is expected at the thermal temperature 
of the flares. However it is possible that the source of hard radiation con- 
tains a hybrid of thermal/non thermal electrons at which point a break is no 
longer necessary in the hard spectra. 

2 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/ 

3 http://suzaku.gsfc.nasa.gov/docs/suzaku/analysis/xisresp 



© 0000 RAS, MNRAS 000, 000-000 



4 Reis, Miller, & Fabian 



column density along the line of sight to the source and its known 
latitude. From estimated bolometric flux and the assumption that 
the source is radiating at less than 5 percent the Eddington lumi- 
nosity, the authors constrained the mass of J1753 to approximately 
4-16 Mq. This range is in agreement with the value of 12 M@ esti- 
mated by Zurita et al. (2008) based on the orbital period of J1753. 
In the same work Zurita et al. estimated a lower limit for the dis- 
tance of 7.2 kpc. For the purpose of our work we will assume a 
conservative range of 4-16 Mq and 7.2-10 kpc for the mass and 
distance respectively. The inner disk inclination was found from 
the reflection features in the XMM-Newton X-ray spectrum to be 
551? degrees (Reis et al. 2009b). 

Using RXTE and XMM-Newton data, Miller, Homan & Mini- 
utti (2006) showed the presence of a cool (kT ~0.2 keV) accretion 
disk extending close to the ISCO in the low-hard state of J1753. The 
presence of this cool accretion disk was later confirmed with RXTE 
observations during its 2005 outburst (Ramadevi & Seetha 2007). 
Reis et al. (2009b) have recently re-analysed the XMM-Newton and 
RXTE spectra of J 1753 in the low-hard state and found that the disk 
likely extends to n»= 3.lio'gr g . This strong constraint was based 
on the self-consistent modelling of both the reflection features and 
the soft thermal-disk component. In this work we further explore 
the spectra described in Reis et al. (2009b) in the 0.5-10.0 keV en- 
ergy range. 



2.5 GROJ1655-40 

GRO J1655-40 (hereafter J1655) has a compact object with a mass 
constrained to 6.3 ± 0.5 Mq and orbital inclination of 70.2°±1.9° 
(Greene, Bailyn & Orosz 2001). Using radio observations, Hjellm- 
ing & Rupen (1995) found a distance of 3.2 ± 0.2 kpc and an inner 
disk inclination of 85° ±2°. In this work we use an inclination of 
70.2° with a lower and upper limit of 68.3 and 87° respectively. 

Using archival ASCA data of J 1655 in outburst, Miller et al. 
(2005) showed evidence of a highly skewed, relativistic line, and 
suggested an inner accretion disk radius of ~1.4r g , indicative of 
a rapidly spinning black hole. This was later confirmed by Diaz- 
Trigo et al. (2007) using simultaneous XMM-Newton and INTE- 
GRAL observations of J 1655 during the 2005 outburst. Brocksopp 
et al. (2006) followed the spectral evolution during the 2005 out- 
burst using the SWIFT X-iay Telescope (XRT). They note the pres- 
ence of a strong iron line at ~6.4 keV in the single observation of 
the low-hard stahfl however they find no evidence of a thermal 
component. It is interesting that the neutral hydrogen column den- 
sity (JVk) remains at an average of ~6.9x 10 21 cm~ 2 for all obser- 
vations reported by Brocksopp et al. other than the single low-hard 
state observation where the value of JVh drops to ~5.9x 10 1 cm~ 2 
(see their Tables 4 and 5). It is plausible that a lower value of Nn 
can act to mask the presence of a weak soft component when the 
signal to noise level is low (See §[3}. 

Takahashi et al. (2008) have recently reported on the Suzaku 
observation of J 1655 in the low-hard state where they find strong 
evidence for both a broad iron line and a soft disk component with a 
temperature of «0.2 keV. Suzaku observed J 1655 in 2005 Septem- 
ber 22, 07:32 UT for approximately 35 ks (Obs ID 100029010). 
The XISs were operated in the 3x3 and 2x2 editing mode with the 



4 Brocksopp et al. (2006) analysed twenty SWIFT pointing covering the 
evolution from the low-hard, high-soft and very-high state, however only 
one pointing captured the source in the low-hard state. 



"1/8 window" employed. XIS0 was not operational during the ob- 
servation. The data were reduced in the standard way as described 
above for CygX-1. Source events were extracted from a circular 
region 70 pixels in radius centred on the source. The background 
is taken from a source-free region 50 pixels in radius towards the 
edge of the chip. The HXD was operated in its normal mode. The 
appropriate response and background files for XlS-nominal point- 
ing were downloaded and the HXD/PIN data were reprocessed in 
accordance with the SUZAKU Data Reduction Guide. The XIS and 
PIN data were fit in the 0.7-9.0 and 12.0-70.0 keV energy range re- 
spectively. A normalisation constant of 1.16 was applied to the PIN 
data in accordance with the latest Suzaku calibration report. 

2.6 XTEJ1118+480 

Dynamical measurements of XTE Jl 118+480 (hereafter J1118) 
have set a strong constraint on the mass-function of 6.1 ± 0.3 Mq 
(Wagner et al. 2001; McClintock et al. 2001). Frontera et al. (2001) 
suggests that the black hole in J 1 1 1 8 has a mass of 7-10 Mq 
with an orbital inclination in the range of 60-80 degrees. A recent 
study by Gelino et al. (2006) places a more stringent constraint of 
8.53 ± 0.60 Mq and 68° ±2° for the mass and inclination respec- 
tively. The authors find a distance of 1.72 ±0.1 kpc to the system. 
To our knowledge the broadest constraint on the distance to Jl 1 18 
was that imposed by McClintock et al. (2001) of 1.8 ± 0.6 kpc. 
Since this range encompass all the latest predictions we employ 
this value in the work that follows. The inclination and mass is 
hereafter assumed to lie in the range 60-83 degree^ and 7-10 Mq 
respectively. 

Jl 118 was observed in its LHS by Chandra, RXTE, UKIRT, 
EUVE, and the HST in 2000 as part of a multiwavelength, multi- 
epoch observing campaign. Based on these observation, McClin- 
tock et al. (2001) reported an apparent cool thermal component 
at «24eV which was interpreted as originating in an accretion 
disk truncated at a radius greater than 70r g . Reis, Miller & Fabian 
(2009) have recently re-analysed the Chandra and RXTE data and 
shown that the accretion disk is consistent with extending close 
to the ISCO, contrary to previous results. The Chandra and RXTE 
spectrum of Jl 1 18 in the low-hard state, as described in Reis, Miller 
& Fabian (2009) is used in the analyses that follows. 

2.7 IGRJ17497-2821 

Unfortunately very little is known about the properties of the newly 
discovered X-ray binary IGR J17497-2821 (hereafter J17497). 
Paizis et al. (2007) have placed the source in the direction of the 
Galactic centre and argue that it must be either very close to or be- 
yond the Galactic centre. We assume the distance to J 17497 to lie in 
the range 5-10 kpc. We are not aware of any estimates for the mass 
or inclination which we assume to be in the range 5-20 Mq and 
10-80 degrees respectively. Paizis et al. (2008) presented the broad- 
band Suzaku observation of J17497 in its low-hard state. They find 
that the data is best modelled with the sum of a thermal disk com- 
ponent plus Comptonisation. The accretion disk is found to have 
an inner temperature of ~0.25 keV and a radius of ~ 123km for an 
assumed distance of 8 kpc and inclination of 60 degrees. For a 10 
solar mass black hole the authors estimate that the disk extends to 
~ 1 1 r g after correction for Compton up-scattering. 

5 http://www.astro.isas.ac.jp/suzaku/analysis/hxd/ 

6 The high value reflects the upper limit predicted by Wagner et al. (2001). 
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Figure 1. Data/model ratio to an absorbed powerlaw. A simple powerlaw does not provide a satisfactory model to any of the sources above. The excess and 
curvature below fa 2 keV is indicative of a thermal disk component which is required at over 5cr level in all sources (see §3.1). GX339-4: EPIC-MOSi and 
2 data are shown in black and red respectively. Cygnus X-l: Obs 1 and 2 are shown in black and red respectively. J1655-40 and J17497-2821: FI and BI 
spectrum are shown in black and red respectively. Jl 118+480: Plus and minus first order spectrum are shown in black and red respectively. J1817-330: Obs 1 
and 2 are shown in black and red respectively. All spectra have been re-binned for plotting purposes only. 
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Table 1. Summary of the various black hole binaries physical parameters. 



Source 


Inclination 


Distance 


Mass 


Total flux 


Disk flux" b 




\ UCglCCft J 


\ K PM 


\ ivl ) 




2 20+° °! 

o.08 


GX 339-4 


10-30 


6-10 


10-20 


16.6 ± 0.1 


J1650-500 


47-70 


1.9-3.3 


5.3-11.3 


213.0 ±2 


61.2 ± 1.2 


Cygnus X-l (1) 


25-50 


2.0-2.2 


7-25 


133.1±f 7 5 


no 7+3.2 
zo - ' -3.0 


Cygnus X-l (2) 


25-50 


2.0-2.2 


7-25 


109.6t|;£ 


15 5+ 3 ' 2 


J1753.5-0127 


49-57 


7.2-10.0 


4-16 


4.10 ±0.01 


0.13 ±0.01 


J1655-40 


68.3-87 


3.0-3.4 


5.8-6.8 


9.5 ±0.2 


0.8 ±0.2 


Jl 1184480 


60-83 


1.2-2.4 


7-10 


13.0 ±0.1 


1.11 ±0.05 


J17497-2821 


10-80 


5-10 


5-20 


15.611:1 


6 2+ 1 - 7 


J1817-330 (1) 


10-80 


1-15 


4-15 


7.15 ±0.15 


2 Q5+ - 25 


J1817-330 (2) 


10-80 


1-15 


4-15 


1 16+ 008 


0.19 ±0.06 



Note - Physical parameters for the various black hole binaries treated in this work. For references and explanations on the range quoted above please see ^2] 
a Unabsorbed flux in units of X 10~ 10 erg cm -2 s" 1 in the 0.5-10 kcV energy range obtained using the XSPEC model CFLUX convolved with the 
model shown in Table 3. b The quoted errors for the fluxes are at the 90 per cent confidence. 



Table 2. Results of fits to simple absorbed powerlaw model. 



Source 


Nn (xlO 22 cm" 2 ) 


r 




N PL C 


x> 


GX 339-4 


0.420 ±0.002 


1.800 ±0, 


004 


0.256 ± 0.001 


6723.5/1635 


J 1650-500 


0.549 ± 0.002 


2.562 ±0 


006 


6.12 ±0.03 


7937.3/1275 


Cygnus X-l (1) 


0.339 ± 0.003 


1.685 ±0, 


005 


1.66 ±0.01 


1088.9/724 


Cygnus X-l (2) 


0.346 ± 0.004 


1.669 ±0. 


006 


1.41 ±0.01 


831.7/685 


J1753.5-0127 


0.175 ± 0.001 


1.666 ±0 


003 


0.0617 ±0.0002 


3044.6/1499 


J1655-40 


0.525 ± 0.003 


1.660 ± 0, 


005 


0.129 ±0.001 


1733.1/1441 


Jl 118+480° 


0.0080+ 0001 


1.877 ±0. 


005 


0.222 ±0.001 


5451.8/4248 


J17497-2821 


4.349 ± 0.003 


1.50 ±0. 


01 


0.114 ±0.002 


1246.6/1184 


J1817-330 6 (1) 


0.12(/) 


3.01 ±0. 


01 


0.215 ±0.004 


403.6/209 


J1817-330 6 (2) 


0.12(/) 


2.01 ±0. 


07 


0.0206 ± 0.0001 


107.7/81 



Note - a A lower limit on the absorbing column density of 8 X 10 19 cm -2 is imposed. fc For Jl 8 17 Nn is frozen at the value indicated. c The powerlaw 
normalisation is referred as N and is given in units of photons keV" 1 cm -2 s~ 1 at 1 keV. All errors refers to the 90 per cent confidence range for one 

parameter of interest. 



We analyse the 2006 September 25 observation (Obs ID 
901003010). The XIS was operated in the 2x2 and 3x3 editing 
mode with the "1/4" window employed due to the brightness of 
the source. The data were reduced in the standard way as described 
above for Cygnus X- 1 . Source events were extracted from a circular 
region 1 10 pixels in radius centred on the source. The background 
is taken from a source free region 110 pixels in radius towards the 
edge of the chip. HXD data is not used in this analysis. Following 
Paizis et al. (2008) we do not use XIS data below 1 keV and limit 
the upper energy to 9 keV. 



2.8 J1817-330 

We are not aware of any constraints on the physical parameters 
of the newly discovered X-ray binary XTE JT817-330 (hereafter 
J 1 8 1 7; Remillard et al. 2006). For this reason we assume the dis- 
tance, mass and inclination to be in the range 1-15 kpc, 4-15 Mq 
and 10-80 degrees respectively. 

The XMM-Newton and INTEGRAL observations of the source 
in the soft state have been reported by Sala et al. (2007). The au- 
thors find that the data can be successfully modelled with the pres- 
ence of a disk component with a temperature varying from approx- 
imately 1 .0 keV to 0.7 keV and an inner radius of ~ 50 km, assum- 
ing a distance of 10 kpc and 45 degrees inclination. The source was 
later tracked through the decline of its outburst by the SWIFT satel- 



lite. Rykoff et al. (2007; hereafter R07) reported on the twenty-one 
snapshot observations tracking J1817 from the high-soft to the low- 
hard state. They found that the source follows closely the Lx oc T 4 
relation during the decline and argued that the accretion disk does 
not recede after the state transition down to accretion rates as low 
as O.OOlZ/Edd- This work was later challenged by Gierlinski et al. 
(2008) who, using the same SWIFT data as R07, argued that ir- 
radiation into the disk gives radii which are consistent with disk 
truncation. 

Here, we re-analyse the data presented in R07 using the same 
reduction procedures detailed by the authors. Out of the twenty- 
two observations (see Table 1 of R07) thirteen are in the soft state, 
two are in the intermediate state and the latter six observations are 
in the low-hard state as described in Gierlinski et al. (2008). In 
this work we will use observations 17-18 since the data quality 
in the latter observations of the low-hard state does not strongly 
require the presence of a disk. Following R07, observations 17 
and 18 (hereafter Obs. 1 and 2) were fit in the 0.6-10.0 and 0.3- 
10.0 keV respectively with the absorbing column density fixed at 
1.2xl0 21 cm" 2 . We note that the neutral hydrogen column den- 
sity in the line of sight to the source is not thought to vary between 
observations (Miller, Cackett & Reis 2009). Table 1 summarises 
the range in the various physical parameters used throughout this 
work. 
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Table 3. Results of fits with the MCD component DISKBB and a powerlaw. 



Source 
GX 339-4 
J 1650-500 
Cygnus X-l (1) 
Cygnus X-l (2) 
J1753.5-0127 
J1655-40 
Jl 118+480 
J17497-2821 
J1817-330 (1) 
J1817-330 (2) 



iV H (xlO 22 cm~< 
0.495 ± 0.006 
0.556 ± 0.004 
0.53 ± 0.02 
0.50 ± 0.02 
0.197 ±0.004 

0.63 ± 0.02 
0.022 ± 0.003 
4.72 ± 0.08 
0.12(f) 
0.12(f) 



1.665 ± 0.009 
2.10 ± 0.01 
1.71 ± 0.01 
1.70 ±0.01 
1.61 ± 0.01 
1.67 ±0.01 
1.69 ± 0.01 
1.56 ± 0.01 
2.1 ± 0.1 
1.5 ±0.2 



n 999 +0.003 
u.zzz_ 002g 

3.40 ± 0.06 
1.752 ±0.025 
1.485 ±0.026 
0.0571 ± 0.0006 
0.1333 ±0.0019 
0.1864 ±0.0028 

n 19H9+0 0034 
U.12»2_ 0032 

097+ - 013 
UiUy( -0.012 

0.012 ±0.003 



kT{keV) 

0.254 ± 0.006 

0.310 ± 0.004 
104+0.005 

u - 13i -0 .004 
104+0.007 

274+ ' 015 
u - z '^-0.014 

0.21 ±0.01 
0.21 ±0.01 
0.20 ±0.01 
0.20 ±0.01 
0.21 ±0.04 



biskbb xlO 3 

5 02+ ' 80 
°' uz -0.67 

55 ±4 
236±|| 

155^^g 

32+ ' 11 
U - J -0.08 

5 4+ 2 - 7 
°-*-2.0 

'•*-1.2 

27±l 

1 3+ 15 
-■-• J -o.6 



Inner radius bc 

2 7+ 2 ' 1 

7 3+ 7 - 1 
'•°-3.3 
5 fi+7-0 

4. 6 +6.3 



-2.7 

^+2.3 
-0.9 

5.9±S:1 
9 Q+2.6 

,1+40 
-12 
1+49 
-13 
3+12 
-3 



1.5 1 



14" 1 
1^ 
3^ 



2874.0/1633 
1507.9/1273 
783.1/722 
719.2/683 
1961.0/1497 
1618.8/1439 
3747.3/4246 
1102.5/1182 
204.0/207 
69.1/79 



Note - Powerlaw nomalisation in units of photons keV - 1 cm -2 s _1 at 1 keV. b The errors quoted for the various inner radii are the 3cr errors estimated 
using Monte Carlo simulations (see §3.2). c The inner radii are given in units of r g = GM/c 2 . 
All errors refers to the 90 per cent confidence range for one parameter of interest. 

Table 4. Same as Table 3 but for the model containing EZDISKBB plus a powerlaw 



Source 
GX 339-4 
J1650-500 
Cygnus X-l (1) 
Cygnus X-l (2) 
J1753.5-0127 
J1655-40 
Jl 118+480 
J17497-2821 
J1817-330 (1) 
J1817-330 (2) 



0.495 1 



-+0.006 
-0.005 

0.556 ± 0.004 
0.52 ± 0.02 
0.50 ±0.02 
0.197 ±0.004 

0.63 ± 0.02 
0.022 ± 0.003 
4.72 ± 0.08 
0.12(f) 
0.12(f) 



1 666+ 009 
2.10 ±0.01 
1.71 ±0.01 
1.70 ±0.01 
1.609 ±0.007 
1.67 ±0.01 
1.69 ±0.01 
1.56 ±0.01 
2.1 ±0.1 
1.5 ±0.2 



0.222 1 



,+0.003 
-0.002 

3.42 ± 0.06 
1.751 ± 0.025 
1.484 ± 0.026 
0.0572 ± 0.0006 
0.133 ± 0.002 
0.186 ± 0.003 
0.128 ± 0.003 
0.0981°;™ 
0.012 ± 0.003 



kT(keV) 

0.241 ±0.005 

0.293 ±0.004 
n iot;+0-005 

U-l»O_0.004 

n iaK4_+0-007 

U.lS5±_u.0 06 

0.26 ±0.01 

0.20 ±0.01 

0.20 ±0.01 

0.19 ±0.02 

0.19 ±0.01 

20+ - 04 
u ' zu -0.03 



^Ezdiskbb X 10 3 
l.lj 
12.0^ 
49j 
32j 
0.075 ± 0.02 



,+0.2 
-0.1 

-1+0.9 
-0.8 
-13 
-11 
-13 
-10 



1 1+O.6 
1-1 -0.4 



1.6 

10.7 
5.9 
0.27 



+0.3 
0.2 
+9.4 
4.7 
+1.9 
1.3 

+0.32 



Inner radius 

q cr+4.4 
O.J_ 2 2 

g+14 
9 -6 
7+13 
'-5 

7A +12.0 



-0.13 



3.0 

17 
18 
1 



4.9 
+4.6 

1.0 
+57 
15 
+72 
17 
+ 17 
3 



X 2 /" 
2885.5/1633 
1509.1/1273 
782.4/722 
718.9/683 
1969.2/1497 
1618.9/1439 
3746.0/4246 
1102.4/1182 
204.5/207 
69.1/79 



3 ANALYSIS AND RESULTS 

The X-ray spectrum (0.1-100 keV) of stellar mass black hole bina- 
ries in the low-hard state above a few keV is usually thought to be 
dominated by inverse Compton scattering of the soft thermal-disc 
photons in a cloud of hot electrons or "corona" (Zdziarski & Gier- 
linski 2004). The geometry of the Comptonisation region is still not 
firmly established. Certain models have hot electrons confined in 
magnetic flares above the cool accretion disk (Beloborodov 1999; 
Merloni, Di Matteo & Fabian 2000; Merloni & Fabian 2001) whilst 
others predict that the hard X-ray emission originates in the base of 
a jet (Markoff & Nowak 2004; Markoff, Nowak & Wilms 2005). A 
further possibility is that the accretion flow consists of a thin disk 
truncated at large distances from the black hole (R tr > 10 3 r g ; 
Esin et al. 1997) and the hard emission originates in the hot, quasi- 
spherical advection-dominated accretion region close to the central 
source. In all these cases the X-ray continuum in the ~2-10keV 
energy range can be phenomenologically modelled by a simple 
powerlaw with a photon index between approximately 1 .4-2, with 
major differences between the models arising at energies greater 
than ~30 and below ~2 keV. 

We begin our analyses by fitting the various spectra up to a 
maximum energy of lOkeV with a powerlaw modified by absorp- 
tion in the interstellar medium (PHABS model in XSPEC). The 
equivalent neutral hydrogen column density (Ng) is initially a free 
parameter in all observations. It is clear from Table 2 and FigQ] 
that the model is a poor description of the data to the majority of 
the objects. Strong curvature below ~2keV is clearly seen in all 



10 r 



£ 

8 
2 



Pi 



1.3 







= A 1 1 







Energy (keV) 



10 



Figure 2. Ratio to powerlaw for Cyg X- 1 when the column density ivjj is 
constrained to 5-8x 10 21 cm~ 2 . See 53. 



sources other than J17497 and J 18 1 7 possibly due to the poorer 
low-energy coverage in those observations. The excess and curva- 
ture seen at soft energies is characteristic of a thermal disc emis- 
sion and cannot be properly modelled by variation in the absorbing 
column density. This is a consequence of the high quality data pre- 
sented here. As an example of the way artificial variation in the 
absorbing column density can act to mask out the presence of a 
disk component we constrain the value of JVh for Cyg X-l to 5- 
9 xlO 21 cm -2 ; this covers the largest published range of TVh for 
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this object (Dotani et al. 1997; Takahashi et al. 2001; Miller et al. 
2002; Torres et al. 2005; Makishima et al. 2008). Fitting the spec- 
tra with a powerlaw now results in a much more distinct soft excess 
as can be seen in Fig. [2] A change in Nh from ~3.4xl0 21 cm~ 2 
(Table 2) to 5 x 10 21 cm~ 2 results in a dramatic difference to the 
modelled spectra of CygX-1 in the low-hard state. A similar be- 
haviour is seen in Jl 1 18 where Nh peaks at the lower limit imposed 
of 8 x 10 19 cm -2 . This behaviour has been reported for CygX-1 
by Takahashi et al. (2001) using ASCA data and possibly explains 
the lower value of Nh found in the single SWIFT pointing of J 1655 
in the low-hard state (see Tables 4 and 5 in Brocksopp et al. 2006). 

Evidence of reflection features in the form of an excess at ap- 
proximately 4-7 keV can be seen in all spectra other than J1118, 
J 17497 and J 18 17. In the work that follows we will ignore the 4- 
7keV range in GX 339-4 and J 1650 and the 5-7 keV range in 
CygX-1, J1753 and J1655. The reflection features will be indi- 
vidually explored in Q3.5\ Table 2 details the parameters obtained 
when fitting the various spectra with a simple absorbed powerlaw. 
Although none of the fits are statistically acceptable we show the 
results obtained with this model to emphasise the necessity of a 
further disk component in the spectra of black hole binaries in the 
low-hard state. 



3.1 The requirement for a MCD component 

The most widely used multicolour disk model (MCD) for the study 
of black holes accretion is the DISKBB model (Mitsuda et al. 
1984). This simple model describes the spectrum from a geomet- 
rically thin and optically thick (Shakura & Sunyaev 1973) ac- 
cretion disk consisting of multiple blackbody components. It is 
parametrised by the colour temperature and a normalisation factor 
defined as (r / d) 2 cos6 where r is the inner radius of the accretion 
disk in km, d is the distance to the source in units of 10 kpc and 8 is 
the disk inclination. Table 3 details the parameters obtained when 
DISKBB is used to model the soft disk emission in the various 
sources. The addition of this model is required at greater than the 
5a level in all observations established using the FTEST command 
in X SPEC. The low temperature found in all sources would render 
the strong detection of the disk component impossible with RXTE. 

The multicolour disk model DISKBB assumes a non-zero 
torque boundary condition at the inner edge of the accretion disk 
(Gierlinski et al. 1999). This is contrary to the standard idea where 
the accretion disk extends down to an inner radius r m at which 
point the material free-falls rapidly with negligible viscous interac- 
tion. In order to test the effect of zero-torque and non-zero torque 
boundary conditions, Zimmerman et al. (2005) constructed the 
MCD model EZDISKBB implementing the zero-torque bound- 
ary condition. Similarly to the DISKBB model, EZDISKBB has 
as a parameter the colour temperature as well as a normalisation 
described as (1/ f) 4 (r /d) 2 cos9, where we remind the reader that 
the colour correction factor, / is the ratio between the colour tem- 
perature and the blackbody temperature of the disk, and all other 
symbols are similar to that of DISKBB. 

To confirm that the results presented above do not depend on 
the choice of MCD model used, we have re-analysed the data with 
the model EZDISKBB. In what follows we assume the standard 
value of 1.7 for the colour correction factor (Shimura & Takahara 
1995) and set a lower and upper limit of 1.3 and 2 respectively 
when calculating errors on the derived inner radius (see i]3.4t . Ta- 
ble 4 shows the results of the fits with the DISKBB model replaced 
by the zero-torque model EZDISKBB. Statistically we find no dif- 
ference in the models between most of the fits, with the two models 
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Figure 3. Inner radius versus Total (unabsorbed) luminosity (0.5-10 keV). 
(Top:) r\ n derived from the normalisation of the DISKBB MCD model. 
(Bottom:) Similar to above but for the EZDISKBB MCD model. The hori- 
zontal dashed line shows ri n =\0r g . The errors represent the 3cr confidence 
level estimated using Monte Carlo simulations (see §3.2). 



also giving very similar colour temperatures. We note here that all 
sources in our sample strongly require the presence of a thermal 
component. The question we will address below is whether this 
thermal emission arises from a small region extending close to the 
ISCO or from a larger, cooler region truncated far from the black 
hole. 



3.2 Inner radius from the thermal component 

The observed thermal disk flux depends on the emitting 
area (oc n n ), inclination and distance to the source. These 
are parametrised in the normalisation of the DISKBB and 
EZDISKBB component as described above. Using the constraints 
on the physical parameters of the various sources along with the 
measured normalisations we derived constraints on the inner-most 
radius of the disk for the various sources. The derived values for 
the inner-radius of the various sources in units of gravitational radii 
are shown in Tables 3 and 4. 

Figure 3 shows the derived inner radii from the normalisations 
of the DISKBB (Top) and EZDISKBB (Bottom) model as afunc- 
tion of the total unabsorbed luminosity. The inner radii obtained 
from the DISKBB model are shown as a function of disk lumi- 
nosity in Fig. 4. The 3cr errors on the inner radii and luminosities 
shown in Figs. 3-6 were estimated using Monte Carlo simulations. 
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Figure 4. Inner radius versus disk luminosity (0.5-10 keV). The inner ra- 
dius, ri n is derived from the normalisation of the DISKBB MCD model. 
The horizontal dashed line shows r; n =10r g . All errors represent the 3rr 
confidence level. 
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Figure 6. Inner radius obtained from the normalisation of the DISKBB 
model versus disk temperature. The errors on the inner radii and disk tem- 
perature are 3 and 1.64ct respectively. 
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Figure 5. Inner radius versus Lx/LEdd. (Top:) n n from DISKBB, (Bot- 
tom:) r in from EZDISKBB. 



A uniform distribution was assumed for the masses, distances and 
inclinations of the various sources with the range shown in Table 1. 
The normalisation of the MCD components and the fluxes were 
assumed to have a Gaussian distribution. When constraints, how- 
ever mild, on the inclinations, masses and distances are available 
we see from Figs. 3 and 4 that the derived inner radii from the ther- 
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Figure 7. Disk luminosity as a function of disk temperature. Assuming a 
luminosity-temperature relation of i^isfe °c Jgj we would expect to see 
a factor of 2-3 change in temperature for the range in luminosity shown in 
the plot. The errors on the luminosity and disk temperature are 3 and 1.64<r 
respectively. 



mal continuum using the model DISKBB are below ~10r g and 
always consistent with extending to the innermost stable circular 
orbit. In all cases other than the two sources with little constraints 
on the physical parameters (J 17497 and J 1817) the maximum value 
possible for the inner radius lies below ~20r g at the 3a confidence 
level. The results obtained from the EZDISKBB model (Fig. 3; 
Bottom) are similar to that of DISKBB but for a factor of pal.34 
higher. Figs. 5 and 6 show the inner radii as a function of Eddington 
ratio and colour temperature respectively. 

It can be seen from Figs. 3 to 5 that the various black hole 
binaries investigated cover a wide range of both luminosities and 
Eddington ratios, with the least luminous source (Jl 1 18) being ap- 
proximately two orders of magnitude less luminous than J 1650. It 
is interesting to note from Fig. 5 that the majority of sources are 
clustered at ~5-9x W~ 3 LEdd which is indicative of a mild selec- 
tion effect in the sample. 
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Source 
GX 339-4 
J1650-500 
Cygnus X-l (1) 
Cygnus X-l (2) 
J1753.5-0127 
J1655-40 
Jl 118+480 
J17497-2821 
J1817-330 (1) 
J1817-330 (2) 



4Q6+ - 006 
u.^ao_ 005 

0.556 ± 0.004 

0.52 ±0.02 

0.50 ±0.02 

n iqs+ 004 

U-19«_0.005 

0.63 ± 0.02 
0.023 ± 0.003 
4.72 ±0.08 

0.12(f) 

0.12(f) 



1.666 1 



r 

,+0.009 
-0.008 

2.10 ±0.01 
1.71 ± 0.01 
1.70 ±0.01 
1 610+ 006 
1.67 ±0.01 
1.69 ± 0.01 
1.56 ± 0.01 
2.1 ± 0.1 
1.5 ±0.2 



,+0.003 

-0.002 

3.42 ± 0.06 
1.75 ± 0.02 
1.48 ± 0.03 

n rK7-?+ - 0005 
u.uo/o_ 0006 

0.133 ± 0.002 

0.186 ± 0.003 

0.128 ± 0.003 

0.012 ± 0.003 



0.222 1 



kT(keV) 

0.241 ± 0.005 

0.293 ±0.004 
n i«c;+o.oo5 

U-l»O_0.004 
n iocr+0.007 

0.25 ±0.01 
0.20 ±0.01 

0.191 ±0.015 
0.19 ±0.01 



0.20 



h0.04 
-0.03 



N 6 

0.09 ±0.01 

. Qo+0.08 

4 1+ 11 
2 71+ 1 - 1 
0.006 ± 0.002 
096+ - 047 

U.UMO_ Q 03 g 

o is5+o ° 3 

Q+0- 8 

5+ - 2 

02+ - 03 
u - uz _0.01 



X 2 /" 
2887.6/1633 
1508.9/1273 
782.5/722 
718.9/683 
1970.4/1497 
1618.9/1439 
3744.9/4246 
1102.4/1182 
204.5/207 
69.0/79 



Source 
GX 339-4 
J1650-500 


N H 
4QO+0.004 

0.557 ±0.004 


r 

1.67 ± 
2.10 ± 






01 
01 


AfpL 
n 999 +0.002 
"'"^-0.003 

3.40 ± 0.06 


kT ( keV) 

254 +0.008 
u.zoi_ 004 

0.311 ± 0.004 


Cygnus X-l (1) 


0.53 ±0.02 


1.71 ± 





01 


1.75 ±0.02 


n iqci+0.005 

U.1SM_q q 04 


Cygnus X-l (2) 


0.50 ±0.02 


1.70 ± 





01 


1.49 ± 0.03 


n 195+0.007 


J1753.5-0127 


0.196 ±0.002 


1.61 ± 





01 


n570 +0 ' 0007 
u - uo ' u -0.0004 


°- 28 -0 1 02 


J1655-40 


0.63 ±0.02 


1.67 ± 





01 


0.133 ±0.002 


0.21 ±0.01 


Jl 118+480 


0.023 ± 0.003 


1.69 ± 





01 


0.186 ±0.003 


0.210 ± 0.007 


J17497-2821 


4.72 ±0.08 


1.56 ± 





01 


0.128 ±0.003 


0.20 ±0.02 


J1817-330 (1) 


0.12(f) 


2.1 ± 





1 


097+0- 013 

u - uy( -0.012 


0.20 ±0.01 


J1817-330 (2) 


0.12(f) 


1.5 ± 





.2 


0.012 ±0.003 


0.21 ±0.04 



N 100 
0021 5+0-0002 

u.uuzio_ 0004 
0.022 ± 0.002 



in+0- 03 

q+0. 03 
-0.02 



0.06" 1 



o oooi 2+ °ooo 5 

u - uuuiz -0. 00002 

0.002 ± 0.001 

0030+ 0006 
u.uuou_ 0005 

02+ - 02 
u - uz -o.oi 

011+ 003 
u - ull -0.002 

0005+ 0006 
u.uuuo_ 0002 



2880.4/1633 
1506.7/1273 
783.42/722 

719.5/683 
1963.5/1497 
1618.6/1439 
3742.7/4246 
1102.5/1182 

203.6/207 
69.0/79 



Table 5. Results for DISKPN fits with the inner radius fixed at 6r g (Top) and 100r g (Bottom). It can be seen that for most sources the parameters and statistics 
do not vary between models. The only exception to this are the large differences in the normalisation of DISKPN. We can use this difference in conjunction 
with physical parameters (inclination, mass and distances; see Table[T} to differentiate between these two contradictory interpretations (see §3.4). All errors 
are quoted at the 90 per cent confidence level. 



3.3 Standard disk in the low-hard state 

The standard Shakura-Sunyaev (SS) model predicts the accre- 
tion disk temperature to be ~ (M/10 M ) _1/4 (i/ L E dd) 1/4 keV 
(Shakura & Sunyaev 1973) for a disk extending to the ISCO. This 
gives a range in temperature of ~0. 15-0.35 keV for the observed 
values of L/L E dd (~5 x 1CT 4 to 1.5 x 10~ 2 ; Fig.[5}. It can im- 
mediately be seen from Figs. 6 and 7 that the range in temperature 
observed is in agreement with that predicted by the SS model given 
the values of Eddington ratio found here. This shows that the be- 
haviour of the accretion disk in the low-hard state for the various 
sources observed are characteristic of that predicted for a standard, 
geometrically thin disk. The results presented so far agree with that 
presented by Rykoff et al. (2007) where the authors claim to see the 
presence of a standard accretion disk extending down to the ISCO 
until at "least" W~ 3 L E dd in Jl 8 17. 



3.4 Can a disk at 100r g radii explain the data? 

In the previous section we have shown using two well known mul- 
ticolour disk models that when the physical parameters of the sys- 
tems are used in conjunction with continuum modelling, the accre- 
tion disk likely extends close to the ISCO in all sources investi- 
gated. 

As previously mentioned, the ADAF model for accretion flow 
in the low-hard state has two distinct zones, with the inner part be- 
ing modelled as a hot, optically thin advection dominated region 
while the outer part consists of a standard Shakura-Sunyaev disk. 
The transition radius is usually assumed to be Rt r ~O(10 3 )r g 
(Esin et al. 1997), however Wilms et al. (1999) and Esin et al. 



(2001) have shown that for GX 339-4 and Jl 1 18 this radius is more 
likely to be 200-400r g and «110r g respectively. In the work that 
follows we will systematically investigate the possibility of having 
a disk truncated at 100r g in the sample of systems in the low-hard 
state presented here. 

The model DISKPN (Gierlinski et al. 1999) is a further devel- 
opment of the DISKBB model taking into account the torque-free 
inner boundary conditiorQ. This model has three parameters: The 
maximum colour temperature of the disk (T co ;) in units of keV 
the inner disk radius, r m , and the normalisation which is defined as 
m 2 cos6/(D 2 / 4 ), where D is the distance to the source in kpc and 
all other symbols are similar to that of DISKBB and EZDISKBB. 
Given that the inner radius is now a parameter of the model, we 
perform fits where this value is fixed at both 6 and 100r g . As can 
be seen from Table 5, models with and without disk truncation give 
equally satisfactory fits. The temperature obtained in both cases are 
in remarkable agreement and only the normalisation differ signifi- 
cantly. 

The normalisation of DISKPN, as well as of any other mul- 
ticolour disk model, conveys important physical information and 
can be used to set apart different, often contradictory, interpreta- 
tions. Using the constraints on the masses and distances to the var- 
ious sources (Table 1) together with the DISKPN normalisations 
for both r ln =6 and 100r g we plot the inclination as a function of 
the colour correction factor for each source (Fig. 8). This factor is 
usually assumed to be a constant close to 1.7 for a wide range in lu- 



It was in fact this work that pointed out the non-zero torque nature of the 
DISKBB model. 
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Figure 8. Inclination vs colour correction factor / for the various sources. The red and green curves are for a disk with inner radius of 100 and 6r g respectively. 
The shaded regions shows the 3cr error range. The upper and lower limits on the inclination are given by the solid horizontal lines and the maximum expected 
colour correction factor is show by the solid vertical line. For the model to be physically consistent the shaded region needs to be enclosed within the 
inclination range and to the left of the vertical line. This criteria rules out a disk at 100r g for all sources other than J17497 and J1817. For simplicity only the 
first observations of Cyg X- 1 and J 1 8 1 7 are shown. 
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Table 6. Results of fits with the LAOR line profile and a powerlaw above 3 kcV. GX 339-4 and J1650 were fit in the 3-12 keV range. For J1655 we used 
RXTE up to 20kcV and J1753 covered the full 3-100 kcV range. For Cygnus X-l only the XIS data were used up to lOkcV. The latter fit also included a 
narrow Gaussian line at 6.4 kcV. The upper energy limits were chosen so as to exclude the prominent Compton hump present in most sources. The values of 
A^h are frozen at that shown in Tableland are given in units of 10 21 cm~ 2 . 



Parameter 


N H 


r 


N PL 


E Laor (keV) 


q 


(degrees) 


r\n (r g ) 


N Laor (xl0- 3 ) 


x 2 /u 


J 1650-500 


0.556 


2.17 ±0.01 


n 7f) +0.07 

J -' u -o.oi 


6.97-0.09 


7 O+0.5 
'■°-0.4 


65+1 


1 45+O.O3 


37+2 


1352.2/1405 


GX 339-4 


0.495 


1 671+0 012 


n 9 nc:+0.004 


6 QS+ 04 

D.MO_q 21 


3.2 ±0.5 


10+ 17 


9 4+O.3 
Z - 4 -0.5 


fi +0.3 


1906.4/1786 


Cygnus X-l (1) 


0.53 


1.71 ±0.01 


1.75 ±0.02 


6. 97_o.o4 


4 +0.5 

■^-o^ 


20+ 4 


fi 3+ 14 


6.6 ± 1.0 


647.9/636 


Cygnus X-l (2) 


0.50 


1.70 ±0.01 


-, 40+0.02 


6.97-o.os 


n fi +0.9 
J - D -0.6 


20+ 3 


7 1+2.3 
'• L -1.B 


5.9 ± 1.0 


561.3/549 


J1753.5-0127 


0.197 


1.60 ±0.01 


0.056 ±0.001 


6.4(f) 


9 9+0.4 
z,z -0.3 




z -^-1.2 


0.14 ±0.05 


1450.8/1467 


J 1655-40 


0.63 


1 74+ - 1 


n 144+0.003 

u - i44 -0.005 


6 47+ - 07 
°- 4 '-0.06 


2.7 ±0.3 


86.0 ±0.2 


1 - — 0.1 


9 0+0.2 

^•°-0.4 


921.1/922 



Table 7. Results of fits to the full energy range using the reflection model REFLIONX. Both the Compton hump and the iron line profile are modelled 
self-consistently. The inner radius is found from the blurring of the reflection component. The values of Nh are frozen at that shown in Tableland are given 
in units of 10 21 cm - 2 



Parameter Afjr 

J1650-500 0.556 

GX 339-4 0.495 

J1655-40 0.63 



2.19 ±0.02 

1 7 9 + 0.01 
' -0.02 
1 71+O.OI 
!■ '1-0.02 



£ (erg cm s 1 ) 
820±Jg° 
500l« 
240+^ 



q Tin (.rg) 

7.7±q% 1.3 ±0.1 
3.0 ±0.1 2.1 ±0.3 



2.72 



+0.7 
0.5 



1.38 



-0.6 
-0.1 



8 (degrees) 
70- 
28^ 
87-i 



-1 

3 +2 
-18 



2.5] 



N PL 
-+0.2 
-0.1 

n 1 10+0.001 

U.llS-o.003 



N, 



Reflionx (xlO 
+ 11 
13 

12+ 2 



102Z 



-3 
12l 2 



x 2 /v 

1441.5/1461 
1935.6/1839 
957.7/966 



minosity (Shimura & Takahara 1995). Davis et al. (2005) showed 
that below a colour temperature of ~ 1 keV the colour correction 
factor is indeed a constant close to 1.7. Above this temperature disk 
ionisation leads to a slight increase in / however it is found to be 
consistently below ~ 2.2 (Davis, Done & Blaes 2006; Done & 
Davis 2008). Contrary to this, it was suggested by Merloni, Fabian 
& Ross (2000) that / is a relatively strong function of accretion 
rates varying from Ml. 7— 3. 

The solid red and green curves in (Fig. 8) shows the depen- 
dence of the inclination on the colour correction factor for an ac- 
cretion disk at 100 and 6r g respectively. The errors, shown by the 
shaded regions, are the 3<r errors estimated by Monte Carlo simu- 
lations assuming a uniform distribution of the mass and distance to 
the various sources (see Table 1) and the DISKPN normalisation 
shown in table 5. It can be seen that for most sources a disk trun- 
cated at 100r g requires a high colour correction factor which is, in 
most cases, inconsistent with the upper limit of three (solid vertical 
line) set by Merloni, Fabian & Ross (2000). We show in gl]that 
this value is in fact likely to be below ~ 2.4 and the value of 3 
used here is a conservative upper limit. J17497 and J1817 are the 
only sources where a truncated accretion disk is not ruled out due 
to the weak constraints on their masses, distances and inclinations. 
However, we have seen from the previous section that the data for 
these sources are suggestive of a disk extending to within 10r g . 



3.5 Inner radius from disk reflection 

So far, we have only focused on the thermal component which in 
the low-hard state is predominant below ~2 keV. In the following 
section we will divert our attention to the reflection features present 
in five of the sources in our sample. To this end we will ignore the 
energies below 3 keV in all spectra and, where possible, extend 
the high energy using either RXTE or Suzaku PIN data. A detailed 
analyses of each individual source is presented below. 



3.5.1 J 1650-500 

We have initially added the RXTE PC A data up the 12keV to 
the XMM-Newton spectrum in order to have a clearer view of the 
continuum. The absorbing column density was fixed at the value 
quoted on Table 3. An absorbed powerlaw does not provide a sat- 
isfactory fit (x 2 /^ = 1988.1/1411), with the bulk of the residuals 
coming from the iron line region. 

The standard and most generally used model for a relativis- 
tic line profile around a rotating black hole is the LAOR model 
in XSPEC (Laor 1991). It describes a broad line from an accre- 
tion disc surrounding a rotating Kerr black hole, with an emissivity 
profile described by a power-law of the form e r = r~ q . We use 
this model to account for the broad residuals seen at the iron line 
region. The LAOR model is calculated with the presumption of a 
maximally rotating black hole. Relativistic line models in which 
the black hole spin is a free parameter are now available (KYR- 
LINE and KERRDISK models; Dovciak, Karas & Yaqoob 2004 
and Brenneman & Reynolds 2006 respectively) however we choose 
to use the well established LAOR model as we are not interested 
in the precise value of the inner radius/spin and to allow for a more 
direct comparison with previous work. The outer disk radius in the 
LAOR model was fixed at the maximum allowed value of 400r g 
and the inclination was constrained to the values shown in Table 1 
(47-70 degrees) for J 1650. The line energy was constrained to 6.4— 
6.97 keV. This resulted in a marked improvement over the previ- 
ous fit with Ax 2 /Ai/ = —451.1/5. A better fit is achieved when 
we allow the photon indices of the XMM-Newton and RXTE ob- 
servations to differ by ^ 0.15. Miller (2009) has shown that small 
differences in the value of Y between various X-ray instruments 
are expected and we allow such deviations in all combined RXTE- 
XMM-Newton analyses. The best fit (x 2 /^ = 1352.2/1405) in the 
3-12 keV range is shown in Fig. 9 (Top-left) and detailed in Table 
6. Fig. 9 (Top-right) shows a plot A\ 2 versus inner radius obtained 
with the command STEPPAR in XSPEC. Based on this analyses, 
the accretion disk in J 1650 is found to be within 6r g at more than 
the 5a level of confidence (Fig. 9). 

Figure 9 (Bottom) shows the best fit mentioned above ex- 
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Figure 9. Top-left: Best fit model for J1650-500 in the 3-12 keV energy range (Table 6). The XMM-Newton data are shown in black with RXTE/PCA data 
shown in red. Top-right: Goodness-of-fit variation as a function of inner radius. The horizontal lines mark the 1, 2, and 3cr level. A disk radius greater than 6r g 
is excluded at the 5<r level (not shown). Bottom: Ratio of the best fit model above to the full, extended data. The HXTE (green) flux at 25 keV is normalised 
to match that of the PCA. The data have been re-binned for display purposes only. 




tended to 100 keV after the addition of HXTE data. The power- 
law indices of both RXTE instruments are tied to each other. To 
make this figure we have normalised the HXTE data so that the 
flux at 25 keV matches that of the PCA. The large excess peak- 
ing at ~30 keV is the "Compton hump" associated with reflection 
of hard X-rays by a cool accretion disk. In order to model the 
full spectra we will use the reflection model REFLIONX (Ross 
& Fabian 2005) blurred with the same kernel used in the LAOR 
model (KDBLUR in XSPEC). The best REFLIONX fit cover- 
ing the full energy band is shown in Fig. 10 (Left) and detailed in 



Table 7. When using the reflection model REFLIONX the Fe-Ka 
line, as well as the Compton hump are modelled self-consistently. 
It should be noted that the refection component cannot account 
for the soft-excess explored in the previous chapters in any of the 
sources investigated here. To emphasise this point we show in Fig. 
10 (Right) the data/model ratio of the best fit reflection model of 
J 1650 extrapolated to the soft energy range. It is clear from this 
figure that a soft disk component is still required. 
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Figure 11. Top-left: Best fit model for GX 339-4 in the 3-12 keV energy range (Table 6). The XMM-Newton and RXTE-PCA data are shown in black and red 
respectively. Top-right: Data/model ratio for the above model extended to 100 keV. The HXTE (green) flux at 25 keV is normalised to match that of the PCA 
(red). Bottom: Best fit reflection model (Table 7). 
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Figure 12. Goodness-of-fit variation as a function of inner radius for 
GX 339-4. The horizontal lines mark the 1, 2, and 3<r level. A disk radius 
greater than 6r g is excluded at the 5cr level (not shown). 

3.5.2 GX339-4 

The presence of strong iron line emission in the XMM-Newton 
spectrum of GX 339-4 was previously seen in Fig. 1. A power- 
law fit over the 3-12 keV range results in a poor fit with \ 2 l v = 
2617.7/1791. This is dramatically improved by the addition of the 
LAOR line model resulting in Ax 2 /Ap — —711.3/5 with an in- 
ner radius approaching that expected for a maximally rotating black 
hole (Fig. 12). 

Figure 11 (Top-left) shows the best fit model over the 3- 



12 keV energy range. The fit extended to the full 3-100 keV range 
is shown in Fig. 11 (Top-right) with the HXTE flux normalised 
to match that of the PCA at 25 keV. Similarly to the spectra of 
J 1650 this source shows the presence of a strong Compton reflec- 
tion hump as expected from a source with strong fluorescence emis- 
sion. Replacing the LAOR component with the blurred-reflection 
model results in an acceptable fit to the full energy range (Fig. 1 1 
Bottom) with the inner radius again implying a maximally rotat- 
ing black hole. Based on the extent of the gravitational blurring of 
the reflection features in GX 339-4, an accretion disk with an inner 
edge at a distance greater than 6r g is excluded at more than the 5a 
confidence level. 

The value for the inner radius found using the full spectra 
(2.0 ± 0.3r g ) is similar to that obtained in the narrower (3-12 keV) 
range from the iron-line shape alone (2.4Z_o' 5 r g ; Fig 12). This em- 
phasises the robustness of spin measurements made from the shape 
of such strong reflection features. The result presented for the inner- 
most extent of the accretion disk in the low-hard state is similar to 
that found for GX 339-4 in both the intermediate state (2.0-2. 2r g ; 
Miller et al. 2008) and the joint fitting of the low-hard and very- 
high state (2.0-2. 1 r g ; Reis et al. 2008) which confirms that the disk 
remains stable during these three spectral states. 



3.5.3 Cygnus X-l 

The profile around the iron-Ka region in Cyg X- 1 cannot be mod- 
elled with a single relativistic line (Fig 13; Left). However, satis- 
factory fits are obtained for both observations by the addition of 
a narrow Gaussian at 6.4 keV together with the LAOR line. This 
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Figure 13. Left: Data/model ratio for Cyg X-l with a powerlaw and relativistic line. Right With the addition of a narrow Gaussian at 6.4 keV. Obs 1 and 2 are 
shown in black and red respectively. 
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Figure 14. Emissivity versus inner radius contour plot for Cyg X-l. The 68, 90 and 95 per cent confidence range for two parameters of interest are shown in 
black, red and green respectively. It can be seen that for the full range of the emissivity the inner radius is constrained between approximately 5-8r g (Obs 1; 
Left) and 5-9r g (Obs 2; Right) at the 90 per cent confidence level for two parameters . 
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Figure 15. Left: Goodness-of-fit variation as a function of inner radius for Cygnus X-l (Obs 1). Right Same but for Obs. 2. The horizontal lines mark the 1, 
2, and 3<r level. 



resulted in fv = 647.9/636 and ^ jv = 561.3/552 for obser- 
vations 1 and 2 respectively. The Gaussian is required at more than 
the 5o" level of confidence (F-value = 28.0 and Prob. = 1.7 xlO -7 ). 
In both observations the inner radii and the emissivity indices are 
well constrained, and it can be seen from Fig. 14 that the subtle 
correlation between these parameters does not provide a satisfac- 
tory solution with large radii for any physical values of emissivity. 
Table 6 details the various parameters in the model. The inner radii 
are constrained to 6.41^'q and 7.1^.1'^r g (90% confidence) in Obs. 
1 and 2 respectively, as shown in Fig. 15. 



3.5.4 J1655-40 andJ1753. 5-0127 

Figure 16 shows the best fit model consisting of a powerlaw and 
the LAOR line profile for both J 1655 (left) in the 3-20 keV energy 
range and J1753 (right) in the full range. The relatively weak reflec- 
tion features in J1753 does not allow for a strong constraint on the 
inner radius using simple phenomenological models, with a disk 
at greater than 20r g only being ruled out at the 3a level (Fig. 18; 
right). However, Reis et al. (2009b) showed using a self-consistent 
reflection model that when the full energy range (0.5-100.0 keV) is 
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Figure 16. Best-fit model consisting of powerlaw plus relativistic line. Left: J1655-40 in the 3-20 keV range. Right: J1753. 5-0127 in the 3-100 keV range. 
The HXTE data is omitted for display only. 
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Figure 17. Left: Data/model ratio for J1655-40 with a model consisting of powerlaw plus relativistic line fitted in the 3-20 keV range. Right: Best-fit reflection 
model in the full 3-70 keV range 
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Figure 18. Left Goodness-of-fit variation as a function of inner radius for J1655-40. Right Similar to the left but for J1753. 5-0127. 



considered, the value of the inner radius is constrained at approxi- 
mately 2.5-5.0 rg in agreement with the value found in section 3.1 
(Table 3). Extending the fit in J1655 to the full 3.0-70.0 keV range 
again shows the presence of a Compton hump (Fig. 17; left) which 
is successfully modelled with the reflection model REFLIONX 
(Fig. 17; right). The inner accretion disk radius in J1655 is con- 
strained to less than 6r g at more than 5<r confidence (Fig. 18; left). 
Table 7 details the parameters found for the fit over the full energy 
range using REFLIONX. 

In all sources investigated in this section the equivalent widths 
of the Fe-Ka line as modelled with the LAOR line profile have 
been found to be consistently above 70 eV at the 90 per cent confi- 



dence level. The only exception to this was J 1753 where the equiv- 
alent width is found to be approximately 50 eV. In the following 
chapter we consolidate our results and discuss the implication they 
have on models for accretion flow. 



4 SUMMARY AND DISCUSSION 

In this work we have systematically analysed a sample of eight 
high quality, stellar mass black hole spectra obtained with various 
CCD and grating instruments. We show observational evidence for 
the possible presence of an optically-thick accretion disk extending 
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close to the ISCO in all eight sources in the low-hard state. In half 
of the sources the evidence comes from both the highly significant 
thermal component as well as blurred reflection features. 

4.1 Constraints from the thermal disk 

In all the objects studied the presence of a low-energy soft ex- 
cess is seen when the spectra are fit with an absorbed powerlaw 
(Fig.0. This excess is traditionally modelled assuming thermal 
emission from an accretion disk where the flux, and hence nor- 
malisation, is related to the disk inner radius according to N oc 
ri n 2 . In ^3.21 and 13.41 we have used this information along with 
constraints on the inclinations, masses and distances to the sources 
found in the literature to obtain the inner radius of accretion for 
the various sources. The results are shown in Figs.[3}(6] For six out 
of the eight sources investigated the 3a upper limit on their inner 
radii are found to be below ~10r g and in all cases they are con- 
sistent with extending to the inner-most stable circular orbit down 
to luminosities of ~ 5 x 10 -4 L e dd- The lack of constraints in the 
physical parameters (mass, distance and inclination) of J 17497 and 
J 1 8 1 7 results in a large uncertainty in the derived inner radii, how- 
ever even when extreme values are used for these parameters, the 
results are still consistent with lack of disk truncation further than 
~100r g for these sources. All results presented above are robust 
to the choice of multi-colour disk component used. We emphasise 
here that throughout this work we have used the largest range in 
the physical parameters for the various sources obtained from the 
literature. Even with this conservative choice we have shown that 
the 3<r upper limit on the disk radius, as obtained from the flux of 
the thermal component, does not require large disk truncation. 

It is interesting to note that in all sources presented here the 
disk temperature are found to be below ~0.3 keV, with the ma- 
jority clustered at ~0.2keV. The disk flux found at these low 
accretion rates results in an insignificant amount of flux visible 
to RXTE/PCA, with its effective low energy cut-off of ~3 keV. 
Zdziarski et al. (2004) investigated the state transition of GX 339-4 
using RXTE. Their Figure 10 shows a sharp change in the inner ra- 
dius - obtained via the disk flux methodology - during the change 
to the low-hard state. Similar conclusions based on the apparent 
change in the normalisation of the DISKBB model have been 
presented by Belloni et al. (1997) for the black hole binary GRS 
J1915+105 again using RXTE/PCA data above 2.5 keV. 

Using a further multi-colour disk model where the inner radius 
is a free parameter we investigated in § 13.41 the possibility that the 
thermal component originates further than 100r g as expected from 
ADAF models. Figure [8] shows that this conviction was excluded 
for all sources other than J17497 and J1817. Our conclusions de- 
pends on the maximum plausible value for the spectral hardening 
factor /. Hardening arise at higher photon energies where the opac- 
ity of the disk atmosphere is lower and photons emerge from deeper 
layers where gas is hotter. Approximately, / = T/T e g ~ r 1//4 , 
where r is the Thomson depth of the layer from which the photons 
emerge (Ross, Fabian & Mineshige 1992; Davis, Done & Blaes 
2006 in the high r limit). The photons then need to scatter out 
through the outer layers which are cooler than where they formed. 
This introduces a Compton down-scattering break in the spectrum 
at an energy e ~ m c c 2 /r 2 . There are therefore two competing 
effects: seeing to deeper layers making the spectrum harder and 
Compton down-scattering by intervening material making it softer. 
The maximum value of / occurs before the effects balance, which 
is at / = (m c c 2 /e) 1//8 . Therefore if we take e ~ 0.5 keV to mark 
the harder part of the thermal emission, / < 2.4. An absolute limit 



will depend on exactly which part of the spectrum is being fitted. 
We consider that / = 3 is a conservative upper limit. 

The high significance of the thermal component in the various 
sources studied here, coupled with the total flux and temperature 
range (~ 0.15 — 0.3 keV) obtained from multi-colour disk models 
such as DISKBB, strongly argue for the presence of a standard 
optically-thick accretion disk extending close to the ISCO. 



4.2 Constraints from reflection signatures 

Further evidence for an optically-thick disk extending to the ISCO 
in the low-hard state comes from the presence of a broad, skewed 
Fe-Ka line and strong reflection hump observed in five out of the 
eight sources presented. These reflection features arise due to re- 
processing of hard X-ray by the cooler accretion disk (Ross & 
Fabian 1993). In the region surrounding a black hole the strong 
gravitational potential causes these reflection features to become 
highly distorted, with the degree of distortion depending on how 
far from the black hole the emitting region is located. Therefore, 
the shape of the prominent iron-Ka line can give a direct indication 
of the radius of the reflection material from the black hole (Fabian 
etal. 1989, 2000; Laor 1991). 

Section [3~5l shows the reflection signatures of five sources. It is 
clear from Figs 9 to 18 that in all cases investigated the inner radius 
of emission, as obtained from the degree of gravitational blurring 
of these features, is consistent with the ISCO down to luminosities 
of ~ 1.5 x 10" 3 L £ dd- For GX 339-4, J1650, J1655 and J1753 the 
results further imply that the central source is a rotating, Kerr black 
hole with the disk extending to less than 6r g at 5a confidence in 
most cases. This is in agreement with results found from prior stud- 
ies in various states (GX 339-4: Reis et al. 2008; Miller et al. 2008. 
J1650: Miller et al. 2002; Miniutti et al. 2004. J1655: Diaz-Trigo 
et al. 2007; Reis et al. 2009b. J1753: Miller et al. 2006; Reis et al. 
2009b) which supports the idea that the emitting area is not chang- 
ing between these states. A non-rotating, Schwarzschild black hole 
is not ruled out for Cyg X- 1 , however a disk with an edge at greater 
than 20r g is excluded at over 5a. 

For four out of the five sources where an iron-Ka emission 
line is seen, the equivalent width, as measured from the LAOR line 
component, is found to be greater than 70 eV at the 90 per cent level 
of confidence. It was shown by George & Fabian (1991) that for a 
neutral, optically thick accretion disk extending to the ISCO the 
reflection fraction R ~ Wif a /180eV. Beloborodov (1999) later 
showed that reflection fractions as low as ~0.3 can be achieved by 
the mild relativistic bulk motion of a corona away from an accretion 
disk. Based on these arguments it can be seen that Wkc > 60 eV 
is still consistent with originating in an accretion disk extending to 
the ISCO. 



4.3 Radio jets and disk truncation 

The sample investigated here covers a large range of luminosities 
from ~0.05 to 1.5% of the Eddington limit, however this is still 
an order of magnitude larger than the predicted quiescent transition 
value of zz4xW~ 5 L E dd (Gallo et al. 2003) below which the black 
holes are thought to be jet-dominated. The presence of radio-jets 
is usually associated with systems in the low-hard state (see e.g. 
Fender 2001, 2006), and it has at times been attributed to a receding 
(ejected?) inner disk (see e.g. Belloni et al. 1997). In this scenario 
the inner part of the accretion flow in the low-hard state is advection 
dominated. The radio jet is then quenched when the inner region 
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is filled by an accretion disk, as is the case in the high-soft state. 
Our findings suggests that jet production is not initiated at the point 
where the accretion disk starts to recede. Another possibility for 
the production of jets could be due to a change in the ratio of the 
energy dissipated in a corona to that dissipated in the disk between 
the various states, an intrinsic change in the mass accretion rate 
through the disk or even changes in the vertical scale-height close 
to the central black hole. 

It is still possible that below the luminosities studied here, an 
advection flow might be present at which point the disk could be 
truncated. This ADAF solution is usually agreed to be the dominant 
emission process in the Galactic Center source Sgr A* (Narayan, Yi 
& Mahadevan 1995; Narayan et al. 1998) as well as various low- 
luminosity AGNs (LLAGNs; Di Matteo et al. 2000, 2003). How- 
ever, an alternative, jet-dominated accretion flow explanation has 
also been successfully applied to these sources (Sgr A* : Falcke & 
Biermann 1999; Falcke & Markoff 2000; LLAGNs: Falcke et al. 
2000). In order to test whether the accretion disk recedes at lumi- 
nosities below those observed here, as well as the profound connec- 
tion between accretion disk and jets, we strongly encourage deeper 
observations of black holes in the low-hard and quiescent states 
from instruments such as XMM-Newton, Suzaku, and in the future 
1XO. 



5 CONCLUSION 

In this paper we have investigated a sample of stellar mass black 
holes in the canonical low-hard state. By systematically analysing 
their X-ray spectrum we have found that in all cases the accretion 
disk is consistent with being at the inner-most stable circular orbit 
down to luminosities as low as ~ 5 x 10 _4 L_Edd- The main points 
and implication of this paper are summarised as follows. 

(i) In all sources investigated the presence of an accretion disk is 
required at the 5<r level of confidence. The temperature and flux of 
this thermal component is consistent with the L cc T 4 relation and 
with a geometrically thin, optically thick accretion disk extending 
to the ISCO. 

(ii) The presence of reflection features and predominantly an 
iron-Ka emission line with an equivalent width greater than 
~70eV is detected in half of the sample. In all these cases, the 
broadness of the reflection features further suggest that the accre- 
tion disk is not highly truncated. 

(iii) Our findings suggest that transition to the low-hard state are 
driven by changes in the corona (perhaps related to jet formation) 
and not changes in the accretion disk, and 

(iv) Jet production is not initiated at the point where the disk 
recedes. 

(v) Furthermore, we suggest the following strong and weak ob- 
servation criteria for disk truncation: 

• The data must be able to rule-out both a broad iron-Ka line 
with an equivalent width >60 eV AND an effective disk temper- 
ature consistent with L oc T 4 or 

• The data must be able to rule-out either a broad iron-Ka 
line with an equivalent width >60 eV OR an effective disk tem- 
perature consistent with L oc T 4 . 

Whilst the number of sources presented here is small, there is 
a general trend that when the accretion disk is statistically required 
the data suggests that the disk in the low-hard state can remain at 
the ISCO. This result is contrary to the lore that the accretion disk 



in the low-hard state is truncated at hundreds of gravitational radii, 
as required in the strong-advection dominated interpretation. 
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